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ABSTRACT 

We present the results of a survey for intervening HI 21-cm absorbers at intermediate and low 
redshift (0 < 2 < 1.2). For our total sample of 24 systems, we obtained high quality data 
for 17 systems, the other seven being severely affected by radio frequency interference (RFI). 
Five of our targets are low redshift (z < 0.17) optical galaxies with small impact parameters 
(< 20 kpc) toward radio-bright background sources. Two of these were detected in 21-cm 
absorption, showing narrow, high optical depth absorption profiles, the narrowest having a 
velocity dispersion of only 1.5 kms -1 , which puts an upper limit on the kinetic temperature 
of T k < 270 K. Combining our observations with results from the literature, we measure a 
weak anti-correlation between impact parameter and integral optical depth in local (z < 0.5) 
21-cm absorbers. Of eleven Call and Mgii systems searched, two were detected in 21-cm 
absorption, and six were affected by RFI to a level that precludes a detection. For these two 
systems at z ~ 0.6 we measure spin temperatures of T s = (65 ± 17) K and T s > 180 K. A 
subset of our systems were also searched for OH absorption, but no detections were made. 
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1 INTRODUCTION 

Observing atomic hydrogen in absorption is an extremely power¬ 
ful probe of galaxies at all redshifts: the detection sensitivity is not 
dependent on many of the properties of the system under study but 
depends rather on the unrelated characteristics of the background 
source. Detection of the Lyo line with absorption techniques in 
the optical and ultra-violet wavelengths has been very successful in 


the past (e.g., Peroux et al. 

2003 [ Prochaska, Herbert-Fort & Wolfe 

2005, Meiring et al. 2011 

and references therein), providing in- 


formation on the physical state of a wide variety of media: from 
the low density intergalactic medium to high redshift galaxies. The 
damped Ly a (or DLA) systems in particular, with HI column den¬ 
sities in excess of Nm = 2 x 10 20 cm -2 are thought to trace gas¬ 
bearing galaxies and contain most of the neutral hydrogen atoms 
at all redshifts (e.g., |Wolfe, Gawiser & Prochaska||2005| >. Indeed, 
these absorption line systems provide the most detailed informa¬ 
tion on the dense, pre-star-forming interstellar medium in galaxies 
at z > 0. However, DLAs suffer from selection effects associated 
with optically selected background sources against which the red- 
shifted Ly a lines are observed, such as limitations due to the at- 
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mospheric cut-off at z = 1.65 and the potential loss of dusty DLA 
systems due to the selection of quasars in optical colours ( |Vladilo,| 
|Prochaska & Wolfe|20081|Curran et al.|2011| ). 

The HI hyperfine 21-cm line presents several advantages over 
Ly a: it does not saturate and the radio continuum is unaffected by 
dust. Furthermore, the 21-cm line enables a determination of the 
kinematics and gas distribution in the intervening absorbers, which 
should then allow us to directly link these systems to their z = 0 
galaxy counterparts (e.g., |Briggs, de Bruyn & Vermeulen||2001| ). 
Comprehensive reviews of 21-cm absorption line studies of red- 
shifted high column density gas are given by |Kanekar & B riggs 
(2004J and more recently by |Morganti et al. 1(2015) . The combina- 
tion of H121-cm and Ly a absorption spectra yields a measurement 
of the HI spin temperature (e.g., |Wolfe & Davis|1979| >, which is a 
good indication of the kinetic temperature of cold clouds. Hence, 
HI spin temperatures measured over a range of redshifts provide 
information on the evolution of the temperature of dense gas over 
cosmic time (e.g. |Kanekar & Chengalur|2003l |kanekar et al.|2014) . 
There is evidence of a correlation between spin temperature and 
redshift, indicating a lower fraction of the cold neutral medium 
(CNM) in the ISM of galaxies at higher redshifts (z > 2; |Kanekar 
et al.|2014), although |Curran et al.| ( |2005] ), |Curran] ( |2012| > and |Cur- 
ran & Webb ( 2006) have argued that the apparent evolution of spin 
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temperature is largely a geometrical effect. Despite the clear advan¬ 
tages of the 21-cm line, it has one marked complication: there is 
a degeneracy between spin temperature and covering factor, which 
implies that high-resolution radio continuum data are needed to un¬ 
derstand the absorption profile in detail. 

Because of the atmospheric cut-off for Ly a, very few opti¬ 
cally selected DLAs are currently known at redshifts z < 1 (e.g., 
Meiring et al.|201 l||Som et al.|2015| ). 21-cm absorption studies are 
potentially a good means of increasing the sample of low redshift 
DLAs as these studies are relatively successful: there are currently 
~ 20 detections of intervening Hi 21-cm absorption at redshifts 
z < 1 ( [Curran et al.|2010|[Gupta et al.|2013||Ellison et al.|2012[ and 
references therein). A more complete sampling of low-redshift HI 
absorbers will eventually help in accurately determining the evo¬ 
lution of the cosmic mass density of neutral hydrogen 9hi over 
half the age of the universe, especially when the results of large 
scale 21-cm absorption line surveys from the Square Kilometre Ar¬ 
ray and its pathfinders will become available. The star formation 
rate density and the optical luminosity density are much better con¬ 
strained over this redshift range (z < 1), and they are observed 
to decline rapidly (e.g., |Hopkins & Beacom|2006]|Madau & Dick-| 
inson 2014 ). To fully understand the evolution of the stellar and 
interstellar content of the universe, it is essential that 9hi is also 
measured more precisely. Larger samples of HI absorbers also as¬ 
sist a more solid interpretation of the evolution of cosmic metal- 
licity at low redshifts, which is currently based on weight z < 1 
DLAs with good metal abundance measurements |Kulkarni et al.| 
|2005l|Kulkarni et al.|20071|Som et al.|20B) . 

Also on much smaller physical scales low redshift 21-cm ab¬ 
sorption line studies are important. The study of apparent pairs of 
foreground galaxies and background radio-bright quasars can con¬ 
strain the covering factor and clumpiness of cold gas in galaxies 
(see e.g., |Reeves et al.|2015| >. In particular, very long baseline in¬ 
terferometry (YLBI) 21-cm absorption line studies allow measure¬ 
ments of small scale structure in the ISM if the background ra¬ 
dio source is sufficiently resolved on VLBI scales jSrianand et al.| 
|2013| ). In order to find these interesting cases that show spatially 
resolved HI absorption, searches for low-redshift HI absorbers are 
essential. 

In this paper we present the results of a Green Bank Tele- 
scop^J(GBT) survey for new H121-cm absorbers at low and inter¬ 
mediate redshifts. We selected our targets to be either known Ca II 
or Mg II absorbers, or galaxies with small impact parameters to ra¬ 
dio bright background sources. 

A complementary part of our survey consisted of searching 
for OH absorption lines with the GBT at low and intermediate red¬ 
shifts. The systems selected to be observed in the OH 1667 MHz 
line are a subset of those observed in the 21-cm line, concentrating 
on the galaxies with small impact parameters to background radio 
sources. The purpose of this survey was to identify molecular ab¬ 
sorption lines in the centres of galaxies, where it is known that the 
ratio of atomic to molecular gas density is higher (e.g., [Leroy et al.| 
|2009] >. Previous observations of the OH line in known molecular 
absorbers indicate that OH is an excellent tracer of molecular gas 
in intermediate redshift absorbing systems, and its detection also 
enables an accurate H 2 column density measurement ( [Liszt & Lu-| 
|cas|1999[|Kanekar et al.|2005] >. 
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The approach of this paper is to first concentrate on the HI 
part of the survey. Section 2 describes in detail the selection of our 
sample and in section 3 we report on the data taking and data re¬ 
duction. In section 4 we describe the detected HI 21-cm absorption 
systems individually in some detail and in section 5 we discuss the 
statistics of low-redshift HI absorbers. The results and implications 
of the OH absorption part of our survey is described in section 6. 
Finally, section 7 presents the conclusions. Throughout this paper 
we adopt Ho — 75kms _1 Mpc -1 , 9 a — 0.70, and 9 m — 0.30. 


2 HI 21CM TARGET SELECTION 


Three different selection techniques were used to identify candi¬ 
dates for 21-cm absorption observations. We briefly describe all 
three methods here. The complete list of targets is presented in Ta¬ 
ble^ where the column ’selection’ indicates which of the selection 
methods was used to include the target. Note that some targets meet 
the selection criteria of more than one method. 

Our first set of objects consists of Mgll absorbers. It has 
been known for many years that a large fraction of Mg II absorbers 
show 21-cm absorption (e.g. Briggs & Wolfe 1983; Lane, Briggs & 
|Smette|2000| ), and it has been argued that this fraction rises when 
those Mg II absorbers are selected that show the highest equiva¬ 
lent widths. In particular, [Lane, Briggs & Smettej ( |2000| ) claimed 
that almost all Mgll systems with Wo(2796) > 2 A show 21-cm 
absorption and [Curran] |20 1 0| find a weak correlation between the 
21-cm line strength normalised by the column density and the Mg II 
equivalent width. However, Kanekar et al. ( 2009) do not confirm 
a statistically significant relation between Mg II equivalent width 
and 21-cm detectability and report an overall detection rate of 21- 
cm absorption in 25 per cent of the Mgll systems. The lack of 
a strong correlation between Mgll absorption properties and 21- 
cm absorption is somewhat surprising, as it has been shown that 
HI column density measured through Ly a absorption does corre¬ 
late with Mg II equivalent width. [Rao, Tumshek & Nestor] ( |2006| l 
showed that the fraction of DLA systems among Mg II systems 
rises to 42 per cent when those systems are selected for which 
W 0 (2796) > 0.5AandW o (2796)/W o (Fe2600) > 2 andlMenard 


|& Chelouche| |2009| report an 8 a correlation between the mean 
Lya Hi column density and Wo(2796) over the redshift range 
0.5 < z < 1.4. 

In order to increase the sample of Mg II absorbers with H 1 21- 
cm absorption, we cross-correlated Mg II absorbers identified in the 
Sloan Digital Sky Survey (SDSS) with background sources brighter 
than 200 mJy at 1.4 GHz in the FIRST data base ( [Becker, White &[ 
Helfa nd|1995| ). The Mg II sample was drawn from the compilations 
of |Rao, Turnshek & Nestor| ( [2006| ), |York & et al.] ( |2006| >, [Prochter,| 
|Prochaska & Buries] ( |2006| >, and |Sardane, Turnshek & Rao] ( |2014| ). 
Four systems from this list fulfil the DLA criterion of |Rao, Turn-| 
shek & Nestor (2006) and have a background source brighter than 
200 mJy at 1.4 GHz. We added to this sample one absorber from 
the list of |Rao, Turnshek & Nestor] ( |2006| , which is a strong Ly a 
absorber. The total sample of Mg II absorbers is listed in Table [2] 
with the original reference indicated. 

Ca II absorbers form our second group of objects. These sys¬ 
tems probably form a dustier and more chemically evolved subset 
of DLAs. Wild & Hewett < |2005| ) argue that the average E(B — V) 
values and Znll column densities combined with conservative as¬ 
sumptions regarding metallicities and dust-to-gas ratios point to 
Ca II systems having H I column densities above the DLA limit. Se¬ 
lection by Ca II absorption is an effective way of identifying high 
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Table 1. Properties of the targets selected for HI and OH observations. Column 1: the name of the background radio source. Systems for which HI lines are 
detected are indicated in bold face; column 2: J2000 coordinates of the background radio source; column 3: emission redshift of the background source, if 
available; column 4: expected redshift of the absorber. In the case of Mg II, Ca II and Ly a. absorbers, z a b s is the measured absorber redshift and in the case 
of known galaxies the redshift of the galaxy is given. * denotes an associated system; column 5: measured 1.4 Jy flux density as given in FIRST; column 6: 
selection method with impact parameter in kpc g iven if appropriate ; col umn 7: Mgn W2 796 equivalent width in A measured fro m SPSS data; colum n 8: 


similar for Can W q 934 ; column 9: reference (1) 


Chen et al. 


(2)|i 


Churchill 


2001 


, (3) 


Rao, Turnshek & Nestor 


Rao, Turnshek & Nestor 


Rao et al. 


|Sardane, Tur nshek & Rao ( 2014}, (6) |Prochter, Proch aska & Buries ( 2 006) . 

For the Mg II Wq 796 measurements the most recent values are taken from either 
Note that ’small 6’ refers to absorption line systems identified with small impact parameter galaxies, whereas ’SDSS galaxy 
without prior detection of associated absorption lines. 


{2006} , (4) 


York & et al. 


{2006), (5) 


(2006) 


{2011} or|Quideretal.|{201l|. 
refers to low redshift galaxies 


Name 

Coordinates (J2000) 
h m s ° " 

^em 

^abs or z gal 

Si.4 
(Jy) 

Selection 
(impact par [kpc]) 

Mg II W 0 2796 

(A) 

Can W 0 3934 

(A) 

Reference 

FBQS J0154-0007 

01 54 54.36 -00 07 23.3 

1.826940 

1.180 

0.26 

Mg II 

1.399 


6 

SBS 0846+513 

08 49 57.98 +51 08 29.0 

0.583715 

0.0734 

0.34 

SDSS galaxy (19.0) 



this work 

GB6 J0855+5751 

08 55 21.36 +57 51 44.1 


0.0260 

0.62 

SDSS galaxy (9.1) 



this work 

B3 0927+469 

09 30 35.08 +46 44 08.7 

2.032000 

0.6213 

0.31 

Can /Mg ii 

3.255 

0.600 

5 

4C +04.33 

09 46 42.42 +04 19 00.3 


0.0877 

0.55 

SDSS galaxy (20.0) 



this work 

GB6J1103+1114 

11 03 34.79+11 14 43.0 

1.72934 

0.7176 

0.27 

Can 

3.470 

0.307 

5 

Q1148+387 

11 51 29.37 +38 25 52.4 

1.302510 

0.5530 

0.50 

small b/ Mg II (18.9) 

0.482 


2 

B1213+590 

12 16 04.72 +58 43 33.3 

1.453360 

0.7246 

0.38 

Can /Mg ii 

2.786 

0.430 

5 

B1239+606 

12 41 29.59 +60 20 41.3 

2.068660 

1.2381 

0.43 

Can 


1.034 

5 

PG1241+176 

12 44 10.83 +17 21 04.5 

1.273000 

0.5507 

0.43 

small b/ Mg II (20.0) 

0.570 


2 

SBS 1307+562 

13 09 09.75 +55 57 38.2 

1.629370 

0.5096 

0.29 

Mg II 

2.355 


6 

B3 1325+436 

13 27 20.98 +43 26 28.0 

2.084390 

0.9539 

0.66 

Can 

2.020 

0.901 

5 

4C +57.23 

13 54 00.12+56 50 04.7 


0.0955 

0.61 

SDSS galaxy (11.7) 



this work 

J1431+3952 

14 31 20.54 +39 52 41.5 

1.215400 

0.6022 

0.21 

Mg II 

2.474 


6 

SDSS 1445+0347 

14 45 53.47 +03 47 32.5 

1.300080 

1.248 

0.28 

Mg II 

1.01 


6 

PKS 1545+21 

15 47 43.54 +20 52 16.6 

0.264300 

0.2657* 

2.4 

small b/Ly a (11.5) 



1 

PKS 1602-00 

16 04 56.14 -00 19 06.9 

1.628760 

1.325 

1.0 

Mg II 

0.67 


4 

3C336 

16 24 39.09 +23 45 12.2 

0.927398 

0.9310 

2.5 

small b/Ly a (19.0) 



1 

3C336-1 

16 24 39.09 +23 45 12.2 

0.927398 

0.8920 

2.5 

small b/ Mg II (20.0) 

1.622 


2 

PKS 2003-025 

20 06 08.49 -02 23 35.1 

1.457000 

1.2116 

2.1 

Mg II 

2.65 


3 

PKS 2135-14 

21 37 45.17 -14 32 55.8 

0.200300 

0.1996 

3.6 

small b/Ly a (16.5) 



1 

PKS 2149+212 

21 51 45.95 +21 30 13.8 

1.538500 

0.9114 

0.9 

log Ahi =20.7 



3 

PKS 2330+005 

23 33 13.17 +00 49 11.9 

0.169989 

0.1699 

0.35 

SDSS galaxy (17.1) 



this work 

PKS 2355-106 

23 58 10.88 -10 20 08.6 

1.639110 

1.172 

0.77 

Mg II 

1.569 


6 


column densities of neutral hydrogen, and thus samples of DLAs at 
low redshift. The SDSS provides an ideal database of quasar spec¬ 
tra to search for these systems. Here, we use the method described 
by |Zych et al.| < [20(59] ) to select Call absorbers. Five of these sys¬ 
tems are found to be in sight lines towards radio loud background 
sources (> 200 mJy). 

The third group of targets consists of galaxies with small im¬ 
pact parameters to background radio sources. For this list we as¬ 
sembled targets using two methods: i) we assembled galaxies from 
the literature that are identified with known quasar absorption lines 
at intermediate redshift and have a radio loud background source; 
ii) we cross-correlated SDSS galaxies with spectroscopic redshifts 
with the FIRST catalogue to identify quasar galaxy pairs. The im¬ 
pact parameters of the selected systems range from 9 to 20 kpc. 
Table [^summarizes the details of our sample. 


3 OBSERVATIONS AND DATA REDUCTION 

Observations were carried out with the NRAO Green Bank Tele¬ 
scope (GBT) during various observing runs between April 2006 
and September 2007. For the back-end we used the spectral pro¬ 
cessor, with a bandwidth of 1.25 or 2.5 MHz, two polarizations, 
1024 channels and a 5s integration time. This setup results in a 
spectral resolution of approximately 0.5(1 + z)kms -1 or 1.0(1 + 


z)km s 1 . The receivers used for this program were the prime fo¬ 
cus instruments PF1-3 (600 MHz), PF1-4 (800 MHz), and PF2 (1 
GHz), with system temperatures of T sys = 48K, 25K, and 22K, 
respectively, and the L-band receiver with T sys — 20K. The ob¬ 
servations were taken in standard on-off mode, with 5 min on and 
5 min off scans. On-source integration times ranged from 1 to 5 
hours, depending on the flux of the background radio continuum 
source. 

Unfortunately, observations in the frequency range much be¬ 
low 1400 MHz are plagued by intermittent strong radio frequency 
interference (RFI). The targets that could not be observed because 
of strong RFI were usually visited again several days later. For 
four targets we have not been able to acquire any good quality 
data, namely B1239+606, SDSS 1445+0347, PKS 2003-025 and 
PKS 2355-106, at sky frequencies of 634.6, 631.8, 642.2, and 
653.8 MHz, respectively. All these cases of strong RFI occurred 
with observations using the PF1-3 receiver, illustrating the diffi¬ 
culty of achieving high quality data around 600 MHz (correspond¬ 
ing to the HI 21-cm line at z ~ 1.2) with the GBT. 

Data reduction was carried out with the DISH data reduction 
tool, which is based on the aips++ package. DISH allows for a flex¬ 
ible editing of the data in order to remove data sections with RFI 
or obvious correlator problems. Initial data editing was done with 
the autoflag tool, after which all data were scrutinised by eye. Data 
taken in one session were median filtered and a baseline was re- 
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moved by making a third-order polynomial fit to the line- and RFI- 
free channels. Data taken over different days were averaged and any 
residual baselines were removed by means of a final second-order 
polynomial fit. 


4 DETECTED HI LINES 


The strength of an HI absorption feature per observed frequency 
channel can be expressed as the optical depth, defined as 

A S' 

T'm — ln(l- —), (1) 

where S is the background source flux at the observed frequency, 
AS is the absorbed flux, and / is the ’covering factor’, i.e., the flux- 
weighted fraction of the continuum source covered by the absorber. 
In the optically thin case, where AS/S < 0.3, equation [I] reduces 
to 


T ~ 


AS 

fs m 


( 2 ) 


By replacing AS with the rms flux density of the spectrum, 
we can calculate the rms of r, or cr T . We conservatively assume 
that an absorption line is significant if it is detected with a strength 
of 5 a over a velocity width of 10 kms -1 . The integrated optical 
depth corresponding to this detection limit is defined as rdV 5 , 10 . 

The HI column density Nui can be calculated from the ob¬ 
served optical depth using 


N m = 1.823 x 10 18 (T s //) I rdV [cm -2 ], (3) 

where T s is the harmonic mean spin temperature of the absorbing 
gas and / is the covering factor. VLBI continuum observations can 
be used to estimate this covering factor, by calculating the fraction 
of the total measured flux density contributed by the core region. 
The assumption is then that the HI absorption arises against the 

cor^f] 

Eight of our targets have been observed by |Helmboldt et al.| 
( |2007| ) as part of their 5 GHz VLBA Imaging and Polarimetry Sur¬ 
vey (YIPS) of flat radio spectrum sources. Figure [I] shows a repro¬ 
duction of their VLB A maps of these eight targets. The covering 
factors / were estimated from these maps by dividing the flux from 
the central point source component by the total flux measured in the 
map. Since these maps are made at 5 GHz, while our observations 
are at much lower frequencies, we note that the calculation of the 
covering factors is intrinsically uncertain and should be taken as an 
estimate only. The results are tabulated in Table [4] 

Neutral hydrogen 21-cm absorption features were searched for 
by eye in the final spectra. We identify four significant 21-cm ab¬ 
sorbers. In the following we describe these detections is some de¬ 
tail. Table[4]summarizes the characteristics of the reduced HI spec¬ 
tra that do not show any HI absorption. This table lists the values 
of cr r , as well as the upper limits to the integrated optical depth 
rdV 5,10 and the channel width dV in kms -1 . A summary of the 
detections is given in Table [4] 


2 Note that in reality we have no information on the actual extent of the 
absorbing material, and whether it fully covers the background emission 
region. 


Table 2. Covering factors / estimated from 5 GHz VLBA maps from |Helm-| 
|boldt et ah]j2007) presented in Figure]!] 


name 

/ 

SBS 0846+513 

0.97 

GB6 J0855+5751 

- 0.50 

B3 0927+469 

0.55 

B1239+606 

0.66 

PG1241+176 

0.40 

SBS 1307+562 

1.00 

B3 1325+436 

0.99 

J1431+3952 

0.95 


Frequency [MHz] 

875.6 875.8 876 876.2 876.4 



redshift 

Figure 2. Top: GBT spectrum of the Hi 21-cm absorber at z = 0.621 
against B3 0927+469. RFI is severe in this spectrum, as indicated by the 
noisy signatures red-ward of the absorption line. Due to this strong RFI 
we cannot confirm whether the tentative second absorption component at a 
frequency of 875.9 MHz is real. The vertical lines indicate the frequencies 
corresponding to the H I line redshifted to the Mg II and Ca II lines detected 
in this system. The error bars give an indication of the redshift uncertainties 
of these lines; Middle: Optical depth distribution of this spectrum. The peak 
optical depth is 3.2 per cent. A single Gaussian fit to the absorption profile 
is shown as a solid line and the baseline is indicated by a dashed line. The 
scale bar indicates the velocity relative to the redshift corresponding to peak 
optical depth; Bottom: The residual optical depth spectrum after subtracting 
the Gaussian fit. 


4.1 Intermediate redshift HI absorbers 

We identify two significant absorption features at intermediate 
redshifts: the Mg Il-selected absorber at z = 0.602 towards 
J1431+3952 and the Call absorber at z = 0.621 towards B3 
0927+469. The spectra of both detections are presented in Figures 
[2]and[4] where we show the antenna temperature Ta in K in the top 
panel, and the observed optical depth in the lower panel. In Figure[5] 
we show multicolour SDSS images of the two fields, centred on the 
position of the background radio sources. In the next two subsec¬ 
tions the characteristics of these absorbers at z ~ 0.6 are described 
in more detail. 




















HI absorption at low and intermediate redshifts 5 










R.A. Offset (mas) 


Figure 1. VLB A maps of eight of the background sources from our survey. The label ’HI’ in the top right corner of three maps indicates that 21-cm Hi 
absorption has been detected against these targets. The data have been taken from |Helmboldt et aL|(2007) . 


Table 3. Parameters of the spectra in which no 21-cm absorption was detected. 


name 

O’ T 

rdV 5,10 

dV comments 



(km s -1 ) 

(km s -1 ) 


FBQS J0154-0007 

< 0.071 


0.562 

strong RFI near expected line 

SBS 0846+513 

0.0046 

< 0.17 

4.5 

baseline ripple, some RFI 

4C +04.33 

0.0099 

< 0.14 

0.449 


GB6J1103+1114 

0.0069 

< 0.16 

1.75 


Q1148+387 

0.0069 

< 0.14 

0.80 


B1213+590 

0.0049 

< 0.089 

0.89 


B1239+606 



0.57 

lost in RFI 

PG1241+176 

0.0080 

< 0.17 

0.80 


SBS 1307+562 

0.030 


0.78 

RFI exactly at expected line 

B3 1325+436 

0.010 

< 0.21 

1.00 


SDSS 1445+0347 



0.58 

lost in RFI 

PKS 1545+21 

0.0028 

< 0.064 

0.78 


PKS 1602-00 

0.0071 

< 0.20 

2.35 

RFI near expected line 

3C336 

0.0049 

< 0.11 

0.95 


3C336-1 

0.0030 

< 0.068 

1.0 


PKS 2003-025 



0.57 

lost in RFI 

PKS 2135-14 

0.0022 

< 0.052 

0.63 


PKS 2149+212 

0.0021 

< 0.055 

0.99 


PKS 2330+005 

0.012 

< 0.20 

0.75 


PKS 2355-106 



0.56 

lost in RFI but detected by 

Gu 


Note that rdV^io indicates an upper limit to the integral optical depth, assuming a significance of 5cr over a 
velocity width of 10 km s _1 . 


4.1.1 B3 0927+469 

The absorption feature seen in the spectrum of B3 0927+469 has 
a single component with a peak optical depth of 0.032. We can fit 
the absorption line with a single Gaussian with a centre frequency 
of 876.005 MHz, corresponding to a barycentric redshift of z = 
0.6214, and a velocity dispersion of cr = 3.8 ± 0.8 kms -1 . The 
integrated optical depth of the line is rdV = (0.36 ±0.05) km s _1 . 
There is tentative evidence for a second absorption component at a 


frequency of 875.9 MHz, but due to strong RFI we cannot confirm 
this detection. Note that this system was also observed with the 
WSRT by |Gupta et al.[ < |2012| , but not detected, although the upper 
limit on the integral optical depth they quote (rdV = 0.24 km s -1 ) 
is lower than the value we measure. The rms noise in our spectrum 
is approximately 20 per cent lower than what |Gupta et al.||2012) 
achieved. To confirm the reality of our detection, we plot in Fig¬ 
ure [3] separate spectra for each polarisation and for each of the two 
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Table 4. Properties of the detected 21-cm absorbers. 


(J T 


7m ax dV 

(km s _1 ) 


rdV redshift 

(km s -1 ) 


GB6 J0855+5751 
B3 0927+469 
4C +57.23 
J1431+3952 


0.0043 
- 0.0057 
0.0060 
0.0235 


0.24 0.21 

0.032 0.82 

0.14 0.45 

0.22 0.83 


1.02 

0.36 

3.03 

4.29 


0.02581 

0.62146 

0.09525 

0.60185 



Frequency [MHz 


875.6 


875.8 876 876.2 

Frequency [MHz] 


876.4 


Figure 3. GBT spectra of B3 0927+469 for each polarisation product (XX 
and YY), and for each of the observing sessions, separated by three days. 
The absorption feature around 876 MHz is apparent in each of the four 
independent spectra. 


observations of this target, separated by three days. It can be seen 
that the absorption feature is detected in each of the four indepen¬ 
dent spectra. 

This system was first identified as a Ca II and Mg II absorber 
in the SDSS data base. The measured redshifts of these lines are 
z = 0.62126 and z — 0.62139, respectively. The uncertainties on 
these measurements are approximately 30 kms -1 , indicating that 
there is no evidence for a redshift offset between the 21-cm line and 
the Ca II and Mg II lines. From the 5 GHz VLB A maps of |Helm-| 
|boldt et ah] ( |2007| > we estimate that / = 0.55 (see Table [4]). The 
HI column density of the system is therefore 1.1 x 10 18 T S cm -2 , 
which is on the low end of the distribution of known redshifted 
21-cm absorbers ( [Curran & Webb|2006| ). Assuming that Call sys¬ 
tems are DLAs, which implies that Ahi > 2 x 10 20 cm -2 , we 
can put a lower limit to the spin temperature of T s > 180 K. The 
assumption that this absorber has an H I column density above the 
DLA limit is corroborated by the results of |Menard & Chelouche| 
( |2009[ >, who find a strong correlation between the Mgll equiva¬ 
lent width and the H I column density. Using their relation and the 
measurement of Mgll Wq 796 = 3.255 A (Table 1), we find that 
Ahi = (2.4 ± 0.5) x 10 20 cm -2 . Furthermore, Ellison |2006| ) 


and |Ellison, Murphy & Dessauges-Zavadsky| ( |2009| > defined the so- 
called ’D-index’ as the ratio of the Mg II equivalent width and the 
total velocity spread of the MgllA2796 profile xlOOO, to identify 
DLAs. For B3 0927+469 we find D = 8.1, which puts it in the 
zone where the probability of the system being a DLA is 54 per 
cent ( [Ellison, Murphy & Dessauges-Zavadsky|2009| >. 



0.6025 0.602 0.6015 0.601 

redshift 

Figure 4. Top: GBT spectrum of the Hi 21-cm absorber at z = 0.602 
against J1431+3952. The vertical line indicates the frequency correspond¬ 
ing to the H I line redshifted to the Mg II line detected in this system. The 
error bar gives an indication of the redshift uncertainties of this line; Mid¬ 
dle: Optical depth distribution of this spectrum. The peak optical depth is 22 
per cent. A double Gaussian fit to the absorption profile is shown as a solid 
line and the baseline is indicated by a dashed line; Bottom: The residual 
spectrum after subtracting the double Gaussian fit. 
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Figure 5. Multicolour SDSS images of the fields around B3 0927+469 (left) 
and J1431+3952 (right). In both images the background source is in the 
centre and marked with ’Q\ The candidate absorbing galaxies that are dis¬ 
cussed in the text are marked with ’A’. The axes are in arcsec, offset from 
the positions given in Table 1. The colour images are based on a combina¬ 
tion of g—, r—, and i —band images, following [Lupton et al. 1(2004) . 
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The SDSS image of the field around B3 0927+469 in Figure[5] 
shows several possible galaxy candidates that could be responsible 
for the 21-cm absor ption. We used the SDSS photometry a nd the 
Le Phar^] software ( Arnouts et al.|l999 Ilbert et al.|2006| ) to de¬ 
termine photometric redshifts, and found one candidate for which 
the derived redshift is consistent with that of the 21-cm absorption: 
SDSS J093034.22+464426.8 at * = 0.67l°; 2 °. However, the im¬ 
pact parameter would be ~ 150 kpc, which renders it very unlikely 
that this galaxy is hosting the gas that is causing the HI absorption. 


4.1.2 J1431+3952 


After our observations of this system, the same HI absorption line 
was detected in new GBT observations reported by |Ellison et al.| 
The quality of our data and the |Ellison et al.| ( 2012| ) data is 
comparable: they report an rms noise of 4.5 mJy per 0.83 km s -1 
channel, while we measure a noise level of 5.0 mJy. We measure 
a peak optical depth of r = 0.22 in the spectrum of J1431+3952. 
This system shows a double structure, which we fit with param¬ 
eters zi = 0.60189, Z 2 — 0.60181, <j\ = (3.5 d= 0.8) kms -1 , 
a 2 = (4.2 =b 0.8) kms -1 , and peak optical depth values of n — 
0.22, T 2 = 0.19. The peak-to-peak velocity difference between the 
two peaks is 15.3kms -1 . For the integrated optical depth of the 
absorption f eature we measure (4.0 zb 0.3) kms -1 , while 

(2012]) report (3.07 zb 0.34) km s -1 . Also for this system the 


et al. 


redshift of the 21-cm line is consistent with that of z = 0.60215 for 
the Mg II system, given the uncertainty in the redshift of the latter. 

This absorption line system was first identified as a Mg II sys¬ 
tem by |Quider et al.| ( [20TT] ). |Ellison et al.| ( |2012| ) use archival HST 
STIS measurements of the Lya line to measure a column den¬ 
sity of log(iVHi/cm -2 ) = 21.2 zb 0.1. From the 5 GHz VLBA 
maps of |Helmboldt et al.| |2007| we estimate that the covering 
factor is f — 0.95, but the 1.4 GHz maps presented in Ellison 
et al.|(|2012|) s uggest a much lower covering factor of / ~ 0.30. 


Ellison et al.| |2012| claim that the resulting spin temperature of 
T s = (90 zb 23) K is the lowest measured value to date. However, 
using our optical depth measurement, the spin temperature would 
even be lower: T s — (65 zb 17) K. Note that Curran et al. (2007) 


reported a measurement of T s = 60 K in an absorbing system at 
z — 0.656 towards 3C336. 

There is one galaxy candidate to the west of the background 
radio source, at an impact parameter of ~ 5 kpc. Also for this field 
we attempted to use the SDSS photometry to determine photomet¬ 
ric redshifts, but the galaxy candidates are too faint to determine 
any useful redshift estimates. 


4.2 Low redshift HI absorbers 

In addition to the two intermediate redshift absorbers, we identify 
two new strong HI absorbers at very low redshifts. These are a de¬ 
tection at z = 0.026 in the galaxy SDSS J085519.05+575140.7, 
seen against the background source GB6 J0855+5751, and a de¬ 
tection at z = 0.096 in the galaxy SDSS J135400.69+565000.2, 
seen against 4C +57.23. The GBT Hi spectra of these detections 
are shown in Figures [6] and [7] and multicolor SDSS images of the 
two fields are shown in Figure [8] 


Frequency [MHz] 

1384.4 1384.5 1384.6 1384.7 1384.8 1384.9 



redshift 


Figure 6. Top: GBT spectrum of the Hi 21-cm absorber at z = 0.026 
associated with the galaxy SDSS J085519.05+575140.7, seen against the 
background source GB6 J0855+5751. The vertical line indicates the fre¬ 
quency corresponding to the HI line redshifted to the systemic velocity of 
the galaxy. Middle: Optical depth distribution of this spectrum. The peak 
optical depth is approximately 24 per cent. A double Gaussian fit to the 
absorption profile is shown as a solid line and the individual Gaussian com¬ 
ponents are shown by thin solid lines. Bottom: The residual spectrum after 
subtracting the fitted model. 


4.2.1 GB6J0855+5751 

The absorption spectrum against GB6 J0855+5751 shows a very 
narrow profile, which we fit with a Gaussian with central redshift 
z = 0.025810 zb 0.000002, velocity dispersion a = (1.51 zb 
0.09) kms -1 and peak optical depth of r = (0.24 zb 0.01). How¬ 
ever, the profile is asymmetric at the base and a second component 
is required to make a satisfactory fit. We fit this second compo¬ 
nent with an offset velocity of +4.61 kms -1 and a dispersion of 
a = (1.50 zb 0.20) kms -1 . The integrated optical depth in the 
absorption seen against J0855 is rdv — (1.02 ± 0.04) kms -1 . 

The measured width of the absorption feature places an upper 
limit on the amount of thermal broadening of the line. Hence, the 
narrowness of the absorption feature seen against J0855 can be used 
to derive an upper limit on the kinetic temperature of the absorbing 
gas. Turbulent motions could also contribute to the total measured 
width. The upper limit to Tk can be expressed as 

rp ^ ttihAv 2 _ 1.2119 x 10 2 Au 2 
k ^ k B 8 In 2 = 8 In 2 ’ 

where k B is the Boltzmann constant, vtih is the mass of the hydro¬ 
gen atom and Av is the profile width measured at half maximum. 
Using this equation, we find an upper limit of Tk = 275 K for both 
components, which implies that we are probing the cold neutral 
medium (CNM). Now, equipped with a measurement of the inte¬ 
grated optical depth and an upper limit on the kinetic temperature, 
we can use Eq.[3]to place an upper limit on the HI column density 
of 5.0 x 10 20 cm -2 , assuming the covering factor / is equal to 1 
and that the spin temperature is equal to the kinetic temperature. 
The background source GB6 J0855+5751 was tentatively identi¬ 
fied as a Compact Symmetric Object (CSO) using 5 GHz VLBA 


3 


http://lephare.lam.fr 
















































8 M.A. Zwaan et al. 


observations by |Helmboldt et al.| j2007| , which reveal two lobes 
separated by ~55 mas. In a forthcoming paper we present new 21- 
cm spectral line YLBI observations of the target and discuss the 
small scale spatial variations of the absorption. These observations 
show that the absorption arises against both components and the 
integral optical depth varies across the background radio source. 

|Braun||2012| ) recently investigated the correlation between in¬ 
tegrated optical depth and column density from emission and ab¬ 
sorption lines studies from the literature. It was found that the re¬ 
lation between the two quantities, rdv and H I, can be fitted very 
well with a model with a ’’sandwich geometry” of cool gas with 
properties surrounded by layers of warm gas. Interestingly, based 
on our total opacity of rdv = 1.02 this model would predict a total 
HI column density of 6.0 x 10 20 cm -2 , very close to the value that 
we determine above. 

The absorption profile reaches a peak at z = 0.02581, 
which corresponds to a velocity offset of 66 kms -1 blueward 
of the optical redshift (z = 0.02603) of the galaxy SDSS 
J085519.05+575140.7. This galaxy is a blue, low luminosity, low 
surface brightness system (Figure[8]), with an absolute Petrosian r- 
band magnitude of —17.8 mag. The axis ratio as measured by the 
SDSS is6/a = 0.7zb0.1, implying an inclination of i — 45°. The 
angle between the major axis of the galaxy and the vector from the 
centre of the galaxy to the background radio source is 0 = 45°. 
Under the assumption that the 21-cm absorption arises in a rotating 
gas disk, this would imply that the rotational velocity of the disk 
at a radius of 9.5 kpc is equal to 66/ sin(z) cos(0) = 172 kms -1 . 
However, the expected rotational velocity of this galaxy from the 
Tully-Fisher relation is U rot =63 km/s ( [Pizagno et al.|2007] >. There¬ 
fore, we conclude that the cold gas clouds we measure are proba¬ 
bly not part of an ordered rotating gas disk, but rather are part of a 
warped disk or otherwise irregular gaseous structure. 

From the SDSS spectrum of the galaxy, we measure F(Ha) = 
25.3 x 10 -17 ergs -1 cm -2 A -1 , from which we derive a ve ry low 


star formation rate of 0.0025 M©yr 1 , assuming a Chabrier <[2003]) 


initial mass function (IMF) and no correction for dust. We measure 
a metallicity of 12+log(0/H)=8.54, using the N2 index of |Pettini &| 
|Pagel| ( [2004] ). Based on the SDSS photometry, we can estimate the 
stellar mass of this galaxy. We apply the method detailed in|Hatz 


|iminaoglou, Fritz & Jarrett| ( [2009) , which uses the |Bertelli et al. 


( |1994) Simple Stellar Population (SSP) models, and find that SDSS 
J085519.05+575140.7 is a very low massjgalaxy with a stellar mass 
of « 2.0 x 10 9 M© (converted to a Cha brier ( 2003 ) IMF). The stel¬ 


lar mass and metallicity are consistent with the mass-metallicity 
relation of |Tremonti et al.| ( [2004) . 


4.2.2 4C +57.23 

The spectrum against 4C +57.23 has a very different shape, dis¬ 
playing two spectrally separated components. The stronger com¬ 
ponent has a peak optical depth of r = 0.14 and can be fitted 
very well with a Gaussian with a = (5.8 zt 0.3) kms -1 . The 
weaker component reaches a peak optical depth r = 0.038, and 
has a — (9.1 =b 0.5) kms -1 . The integrated optical depth in the 
absorption seen against 4C +57.23 is rdv = 1.17 zb 0.08 kms -1 . 
After subtracting the two components, some residual structure is 
present, in particular in between the peaks of the two components, 
at z = 0.0953. However, the signal to noise in the spectrum is not 
sufficiently high to fit a third component. 

The peak optical depth is seen at a velocity 70 kms -1 
blueward of the optical redshift of the galaxy SDSS 
J 135400.69+565000.2. This appears to be an edge-on spiral 


Frequency [MHz] 



0.0958 0.0956 0.0954 0.0952 0.095 0.0948 

redshift 


Figure 7. Top: GBT spectrum of the Hi 21-cm absorber at z = 0.096 
associated with the galaxy SDSS J135400.69+565000.2, seen against the 
background source 4C +57.23. The vertical line indicates the frequency cor¬ 
responding to the HI line redshifted to the systemic velocity of the galaxy. 
Middle: Optical depth distribution of this spectrum. The peak optical depth 
is approximately 14 per cent. A double Gaussian fit to the absorption pro¬ 
file is shown, and the baseline is indicated by a dashed line. Bottom: The 
residual spectrum after subtracting the fitted model. 
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Figure 8. Multicolour SDSS images of the fields around GB6 J0855+5751 
(left) and 4C +57.23 (right). In both images the background source is in the 
centre and marked with ’Q\ The foreground galaxies can be easily iden¬ 
tified, and are marked with ’G’. The axes are in arcsec, offset from the 
positions given in Table 1. The colour images are based on a combination 
of g—, r —, and i —band images, following [Lupton etaL|(2004) . 


galaxy, probably of morphological type Sa. The SDSS spectrum of 
the galaxy does not show any strong emission lines. The sightline 
to the background radio source intercepts this galaxy along its 
minor axis in a region where some low surface brightness star 
light can be seen. We note that most other galaxies that have been 
seen in 21-cm absorption are of late morphological type, the only 
exception being J0849+5108 at a redshift of z = 0.3, discovered 
by |Gupta et al.| ( [20T3] ). 
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5 LOW REDSHIFT GALAXY-QUASAR PAIRS 

With the two new local 21-cm absorbers reported in this paper, the 
total sample of intervening 21-cm absorption systems arising in low 
z galaxies rises to 14. Table [5] gives an overview of all known low 
redshift systems, listing the redshifts, impact parameter between 
the identified galaxy and the background radio source, and the in¬ 
tegral optical depth rdV. We somewhat arbitrarily adopted a red¬ 
shift limit of z < 0.5 to include the two systems recently pre¬ 
sented in |Gupta et al.|(2013] ), and the system at z = 0.43 identified 
by |Kanekar, Athreya & Chengalur| ( [2502] >. |Gupta et al.| ( [20T3] > re¬ 
ported a weak correlation between impact parameter and rdV, with 
a Spearman rank coefficient of r — —0.3. Adding our two new sys¬ 
tems and two systems from the literature that [Gupta et al.||2013) 
did not consider (J074110.85+311154.7 and KAC2002), we find a 
slightly higher value of r = —0.39, with an associated significance 
level of 0.81. 

Under the assumption that low redshift 21-cm absorbers arise 
in the gas disks of galaxies, such a correlation with a large scatter is 
expected. HI emission line maps of local galaxies typically show a 
declining HI column density toward larger galactocentric radii, but 
for a sample of HI absorbers the distribution of absorption strength 
vs. impact parameter is a convolution of HI profiles and the HI 
size distribution function of galaxies. In |Zwaan et al.|p005) local 
galaxies are used to calculate the probability function of HI cross- 
section in the column density vs. impact parameter plane for ab¬ 
sorption selected systems. The conditional probability distribution 
shows that the median impact parameter b at a column density of 
log Nm = 20.3 (which defines the DLA limit) is 10 kpc, whereas 
the most likely impact parameter is b = 5 kpc. 

It is interesting to compare the sample of low redshift inter¬ 
vening 21-cm absorption systems arising in galaxies with the like¬ 
lihood distribution of HI emission in the Nm — 5-plane. However, 
Equation 2 shows that for the calculation of column densities from 
optical depths, a measurement of the spin temperature is required, 
which is unknown for most of the systems in Table [5] |Chengalur^] 
|Kanekar & Roy|(|2013|) use simul ations to show that the isother¬ 
mal estimator (Dickey & Ben son 1982) can be used to calculate 
HI column densities from optical depth measurements with an ac¬ 
curacy of a factor 2. In order to arrive at a more accurate relation 
between optical depth and column density, |Braun| ( |2012| > used a 
sandwich model of cool gas, surrounded by layers of warm gas. Us¬ 
ing their relation, and the assumption that the total velocity width 
of the absorption profile is 15 kms -1 , we can transform the rdV 
values in Table [5] to column densities. This calculation uses the as¬ 
sumptions that the ’threshold column density’, for which the optical 
depth approaches zero, is No = 1.25 x 10 20 cm -2 and the satu¬ 
ration column density, where the optical depth reaches infinity, is 
Noo = 7.5 x 10 21 cm -2 . 

In Figure [9] we plot the measured impact parameters b against 
these derived column densities. Also shown is the conditional prob¬ 
ability distribution of impact parameter as a function of column 
density as measured by |Zwaan et al.| ( |2005j ) in local galaxies. 
Shown are the median impact parameter and the 25 and 75 per¬ 
centiles. Note that the measured distribution of absorbing systems 
follows the general trend that is expected from a random incidence 
of sight lines through local galaxies. However there is a dearth of 
low impact parameter systems. This effect is expected as quasar- 
galaxy pairs can only be identified if the angular separation be¬ 
tween the galaxy and the background source is sufficiently large 
so that the two sources can be separated. Even if the background 
source is taken from a radio catalogue, the pair would only be iden- 
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Figure 9. Impact parameter vs. derived Hi column density (using the con¬ 
version per [Braun] (2012) ) of low-redshift quasar-galaxy pairs. The ringed 
points indicate the two new low redshift absorbers reported in this paper. 
The lines indicate the probability distribution functions of impact param¬ 
eter derived by |Zwaan et al. | {2005) , based on a large sample of local Hi 
21-cm emission line maps. The solid line is the median impact parameter as 
a function of HI column density, and the long-dashed lines are the 25 and 
75 percentiles. The vertical line indicates the DLA column density limit. 


titled if the separation were larger than several arcsecs, or otherwise 
the background source could possibly be associated with the fore¬ 
ground galaxy. 

|Reeves et al.| ( |2015] ) recently presented a survey for 21-cm ab¬ 
sorption and emission in a sample of six nearby, gas-rich galaxies. 
The impact parameters between the background radio sources and 
the foreground galaxies ranged between 10 and 20 kpc. No HI 21- 
cm absorption was detected in any of their galaxies, while in our 
sample of five galaxies, absorption was detected in two cases. Al¬ 
though the statistics are poor, the difference between our results 
and the |Reeves et al.| (2015) results are probably due to the fact 
that 1) the continuum flux of the background sources in our study 
was higher, and 2 ) our sample includes systems with smaller im¬ 
pact parameters. In fact, we find that out of our sample of five, 
21-cm absorption was detected in the two galaxy-quasar pairs with 
the smallest impact parameters (~10 kpc). |Reeves et al.| ( |2015| ) also 
note that the detection rate of 21-cm absorption is probably strongly 
correlated with the structure of the background source, with com¬ 
pact radio sources producing higher absorption line sensitivity. 

Similar relations between impact parameter and HI column 
density of DLA galaxies have been found at higher redshifts. Using 
optical imaging of a large sample of DLAs, |Rao et al.|(2011) find an 
anti-correlation between impact parameter and HI column density 
at redshifts 0 < z < 1. They find a median impact parameter of 
16.2 kpc for the DLA galaxy sample. [Peroux et aL](201 1) apply the 
integral field unit (IFU) technique to detect redshifted Her emission 
associated with DLAs, and find indications of an anti-correlation at 
redshifts 1 < z < 2. At redshifts z > 2, |Krogager et al.|(2012| ) 
combine their own observations of three DLA systems with liter¬ 
ature results. The combined sample of 10 identifications shows a 
Spearman rank coefficient of r = —0.6 between Nm and impact 












































10 M.A. Zwaan et al. 


Table 5. Summary of all known intervening 21-cm absorption systems arising in z < 0.5 galaxies. 


galaxy 

redshift 

impact parameter 
(arcsec) 

impact parameter 
(kpc) 

rdV 
(km s -1 ) 

reference 

Can W 0 3934 

UGC 7408 

0.0015 

93 

2.6 

1.3 

1 


NGC 3067 

0.0049 

110 

11.1 

0.11 

2,3 

0.43 ± 0.05 

ESO 1327-2041 

0.0180 

38.5 

14 

0.14a 

2,4 

0.58 ± 0.04 

J085519.05+575140.7 

0.0258 

19.6 

9.5 

1.02 

this work 


NGC 6503 

0.0270 

18.7 

10 

0.22 

6,2 

0.35 ± 0.08 

GQ1042+0747 

0.0332 

2.5 

1.7 

0.19 

5 

< 1.07 

0248+430 

0.0520 

15.0 

15 

0.26 

7 

1.52 ±0.17 

J163956.38+112802.1 

0.079 

3.6 

5 

15.7 

8 


J135400.69+565000.2 

0.0952 

6.8 

11 

3.03 

this work 


J124157.26+633237.6 

0.143 

4.4 

11 

2.90 

9 

1.01 ±0.11 

J074110.85+311154.7 

0.2212 

6.0 

20 

0.36 

10,11 


J084958.09+510826.7 

0.3120 

3.0 

13 

0.95 

12 

0.60 ±0.18 

J144304+021419 

0.3714 

< 1 

< 5 

3.4 

12 


KAC2002 

0.4367 

2.2 

11.6 

0.75 

13 



References: 1: 


Borthakur et al. (2 0~j"4); 2:|Carilli & van G orkom (1 992); 3 : |Keeney et al.|(2 005 ); 4: Kee ney e t al. ]M;5: 


Borthakureta .1(2010); 6:|Boisse et al.|(1988); 7:|Hwang & Chiou (2004) ; 8: | Srianand et al. |(20 13 ; 9:|Gupta et al.|(2010) ; 


10 : |Chengalur & Kanekar |(T999) ; ll: |Cohen|(2001) ; 12: |Gupta et al.H2013) 13: |Kanekar, Athreya & Chengalur|(2002y 


parameter, with an associated significance level of 93 per cent (our 
calculation). Using stacking, Noterdae me et al.| <[2014 ) find Ly a 
emission in a sample of extremely strong damped Ly a systems 
with Nm > 0.5 x 10 22 cm 2 at 2 < z < 4. They conclude that 
the impact parameter of these high column density systems must be 
smaller than 2.5 kpc, again confirming the anti-correlation between 
impact parameter and iVm for DLAs. Interestingly, the median im¬ 
pact parameter for DLA galaxies does not seem to evolve strongly. 
While |Zwaan et al.|(|2005|> quote a median impact parameter of 
8 kpc at z — 0, |Krogager et al.| (2012) find a value of 7 kpc at 
z ~ 2.5. Hydrodynamic simulations at z = 3 can reproduce the 
anti-correlation between impact parameter and HI column density 
(Rahmati & Schaye|2014) . 


6 OH ABSORPTION LINE RESULTS 

In addition to our search for H I absorption lines, we observed a 
subset of our targets in the OH 1667 MHz line. For these obser¬ 
vations we selected the group of small impact parameter quasar- 
galaxy pairs and the Ca II absorbers. The motivation to search for 
molecular lines in small impact parameter systems is the observa¬ 
tion that local galaxies usually show CO emission only in their cen¬ 
tral regions. From CO imaging of nearby galaxies (e.g., jHelfer et al.| 
|2003||Leroy et al.|2009| ), we see that the area over which CO is de¬ 
tectable is much smaller than the HI area above the DLA limit. In 
addition, the regions where the CO columns are highest (iVn 2 > a 
few xlO 21 cm -2 ) often show a depression in Hi column density, 
in many cases below the DLA limit. Furthermore, [Liszt & Lucas | 
(1999) demonstrated by comparing Galactic absorption spectra of 
several poly-atomic molecules with that of HI that the molecules 
generally avoid the deepest HI features and often the high column 
density molecules have no corresponding HI. Presumably, in the 
regions where the molecular column densities are highest, most of 
the HI has been converted to H 2 ( |Schaye|2001) . 

|Wild & Hewett| (2005) found that Ca II absorbers have rela¬ 
tively high dust content and typically high HI columns. In a sub¬ 
sequent paper (Wild, Hewett & Pettini|2006) , they showed that the 
Call systems have highly depleted refractory elements and high 
dust-to-metals ratios with values close to, or even larger than, those 


observed locally. These properties suggest that a substantial frac¬ 
tion of the Ca II absorbers are more chemically evolved than typi¬ 
cal DLAs. The combined results of | Wild & H ewett (20 05) ; [Wild, | 
|Hewett & Pettini| (2006) ; [Nestor et al.| ( |2008) ; [Zych et al. ( 2009) 
lead to an interpretation of the Ca II absorbers arising in the inner 
parts of galaxies, preferentially regions of high gas density, with 
cross-sections only 20-30 per cent of those of DLAs. These ar¬ 
guments make them excellent candidates for molecular absorption 
line studies. 

|Liszt & Lucas| (1999) showed that the abundances of OH 
and HCO+ are the best tracers of molecular hydrogen, demon¬ 
strating a near constant relation between their column densities. 
Arguably, OH is a more accurate indicator of H 2 column den¬ 
sity than molecules observable in the millimetre regime, includ¬ 
ing the often used tracer CO. Furthermore, in all five cases where 
molecular absorption has been seen in redshifted absorption line 
systems, corresponding OH absorption has subsequently been de¬ 
tected: B3 1504+377 (z = 0.673) andPKS 1413+135 (z = 0.247, 
Kanekar & Chenga lur||2002) , B 0218+357 (z = 0. 685, [Kanekar 


et al. 


2003), 1830-211 (z = 0.886, |Chengalur, de Bruyn~& 


Narasimha||1999) , and PMN J0134-0931 (z = 0.765, |Kanekar| 


et al.|2005). 


Table [6] lists the systems that were observed in the OH line. 
The observational set-up was the same as that described in section 
3, and data reduction was identical to what was done for the HI 
observations. Unfortunately, none of the systems we observed were 
detected in absorption in the 1667 MHz OH line. We should note 
that seven of the spectra were severely affected by RFI. In particular 
for the low-redshift quasar-galaxy pairs that were observed in the 
frequency range of 1500 to 1600 MHz, no useful data could be 
obtained. 

Table [6] summarises the upper limits to the integral optical 
depth rdV that can be derived from our measurements. To con¬ 
vert the optical depth limits into OH column densities, we use the 
equation given by (Liszt & Lucas|1996) : 


Non = 2.24 x 10 14 x T ex / r 1&m dV [cnT 2 ] 


/’ 


(4) 


where T ex is the excitation temperature in Kelvin. The value of 
T ex is somewhat uncertain. As argued by [Kanekar et al.|(2003| ), the 
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Table 6. Parameters of the spectra in which OH absorption was searched for. No detec¬ 
tions were made. 


Name 

<7 7 - 

TdV 5,10 
(kms -1 ) 

dV 

(kms -1 

comments 

) 

SBS 0846+513 

0.0091 

< 0.20 

0.97 

RFI 

GB6 J0855+5751 




lost in RFI 

B3 0927+469 

0.0045 

< 0.11 

0.71 

RFI 

4C +04.33 




lost in RFI 

GB6J1103+1114 

0.0046 

< 0.11 

1.51 


Q1148+387 




lost in RFI 

8C 1213+590 

0.0032 

< 0.088 

1.52 


B1239+606 

0.0039 

< 0.10 

1.97 


PG 1241+176 

0.0041 

< 0.11 

0.68 

lots of RFI 

B3 1325+436 

0.013 

< 0.33 

1.70 


4C +57.23 

0.0025 

< 0.054 

0.99 


PKS 1545+21 

0.0019 

< 0.051 

1.17 


3C336 

0.011 

< 0.41 

1.66 

strong RFI below 880.6 MHz 

3C336-1 

0.0049 

< 0.13 

1.70 


PKS 213514 

0.00086 

< 0.025 

1.08 


PKS 2330+005 

0.0041 

< 0.081 

1.05 



excitation temperature in redshifted absorbers is typically higher 
than Xcmb (1+z). In the Milky Way, a typical value of T ex = 10 K 
is found in dark clouds. Lacking better constraints, we also assume 
here that T ex = 10 K. Finally, to convert OH column densities 
into H 2 column densities we adopt the relation from|Liszt & Lucas| 
(19991: Nh 2 = 1.0 x 10 * * * * * * 7 x No H . 

In Figure^)] we plot the upper limits to the H 2 column den- 
sides as a function of impact parameter. We also show the mea¬ 
surements in the five Call systems, for which we have no mea¬ 
surement of the impact parameter because the host galaxies of the 
Call absorption have not been identified. Our 1 Vh 2 upper limits 
are typically a few xl0 21 cm -2 . It is interesting to compare these 
values to the cumulative distribution function of impact parame¬ 
ters of H 2 cross-section-selected systems, as derived by jZwaan &] 
Prochaska (2006). There it is shown that, based on the analysis of 
CO maps of local galaxies, it is expected that H 2 column densities 
of log Nh 2 > 21 are only detected if the impact parameter is less 
than 10 kpc. Our non-detections are consistent with these expecta¬ 
tions. 

Also for the Ca II systems at intermediate redshifts, we reach 
upper limits of 1Vh 2 ~ 2 x 10 21 cm -2 . At 2 = 0 we would expect 
to detect these column densities at a median impact parameter of 
3 kpc fZwaan & Prochaska|2006| . As stated earlier, [Wild & Hewett| 

( |2005| interpret the Ca II absorbers as the inner parts of galaxies, 

with cross-sections only 20-30 per cent of those of DLAs. Adopt¬ 

ing a median impact parameter of low-redshift DLAs of 7.6 kpc 
( jZwaan et al.|2005) , we would expect that the median impact pa¬ 

rameter of Call systems is therefore ~ 3.8 kpc. Using this line of 

reasoning, we could expect to have detected OH absorption in one 
or two Ca II systems in our small sample of five. Now that larger 
samples of Ca II absorbers are available, it would be very interest¬ 

ing to repeat the experiment of detecting molecular gas - OH or 
other molecules such as CO - in these systems. 


7 CONCLUSIONS 

We present the combined results of several Green Bank Telescope 
surveys for intervening HI and OH absorption at low and interme¬ 
diate redshifts (0 < z < 1.2). We obtained good quality data for 



impact parameter [kpc] 


Figure 10. Upper limits on the H 2 column densities, derived from the OH 
absorption non-detections. For the local galaxy-quasar pairs, the upper lim¬ 
its are plotted against their impact parameter. The Ca II absorbers are plotted 
at a fixed impact parameter of 4 kpc. 

17 out of 24 systems - observations of the remaining seven sys¬ 
tems were severely affected by radio frequency interference. The 
total sample of observed systems consists of Ca II absorbers, Mg II 
absorbers, apparent quasar-galaxy pairs with small impact param¬ 
eters at low redshift, and Mg II and Ly a absorbers with identified 
host galaxies. Our conclusions are the following: 

(i) Out of five low redshift (z < 0.17) quasar-galaxy pairs, we 
detected strong Hi 21-cm absorption in two systems. The maxi¬ 
mum 21-cm optical depths are 14 and 24 per cent, respectively. 
Both absorption profiles are very narrow, the narrowest has a ve¬ 
locity dispersion of only 1.5 km s -1 , which puts an upper limit on 
the kinetic temperature of T& < 270 K. Combining these two sys- 
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terns with 21-cm absorbers from the literature, we measure a tenta¬ 
tive anti-correlation between impact parameter and integral optical 
depth. Using a sandwich model for the cool gas, we convert optical 
depth measurements into column densities. The resulting correla¬ 
tion between impact parameter and HI column density is in good 
agreement with what is expected if 21-cm absorbers arise in the gas 
disks of galaxies that are seen in 21-cm emission. 

(ii) Of the sample of five Ca II and Mg II systems for which good 
data were obtained, we measure 21-cm absorption in one Mg II sys¬ 
tem, and one Ca II system, both at 2 « 0.6. Since some of the Ca II 
absorbers are also Mg II absorbers, in total three Mg II systems have 
been observed. The detection rate is in agreement with previous 
studies of Mgll systems. Given the indications that Call systems 
are on average higher HI column density absorbers, we would ex¬ 
pect a higher detection rate for these systems, but unfortunately, the 
statistics are too poor to draw any conclusions on this. 

(iii) No OH absorption was detected in any of our systems. 
Comparing our results with the expected strong anti-correlation 
between column density and impact parameter based on what is 
observed in local galaxies, our non-detections are fully consistent 
with the expectations. To increase the sample of known intervening 
molecular absorption lines, even smaller impact parameter systems 
(b < 10 kpc) need to surveyed. When larger samples of Call ab¬ 
sorbers become available, we recommend observations of molecu¬ 
lar absorption lines in these systems with radio or millimetre bright 
background sources. 

By following up known optical quasar absorption line systems 
in the radio, the total sample of 21-cm absorbers has increased 
steadily over the last few decades. However, the next big step in 
the study of 21-cm absorption systems is expected with the ap¬ 
pearance of wide field blind surveys with upcoming instruments 
such as MeerKAT ( [Booth et al.||2009| >, ASKAP ( [DeBoer et al.| 
|2009| ), and APERTIF ( [Oosterloo et al.||2009) . For example, the 
ASKAP FFASlQ survey will search for 21-cm absorption features 
in approximately 150,000 sight-lines to background radio sources 
brighter than 50 mJy at 800 MHz. The full southern sky survey will 
increase the number of known intervening 21-cm absorbers in the 
range Ahi ~ 10 19 cm -2 to 10 22 cm -2 between 0.4 < z < 1.0 
with a factor of 20. The very first pilot results of FFASH ( jAllison| 
|et al.|2015] > resulted in the detection of a new absorber at z = 0.44 
and demonstrate the ability to efficiently perform wide-field HI ab¬ 
sorption surveys. Surveys with APERTIF and MeerKAT will probe 
lower down the column density distribution function, but will cover 
smaller areas. Together, the result of these surveys will lead to a 
much better understanding of the global evolution of cold gas in the 
z <1 Universe, but will also turn up many interesting instances of 
galaxies at intermediate redshifts in which the ISM can be studied 
in great detail. The examples presented in this paper illustrate the 
type of intervening absorption systems that will be detected rou¬ 
tinely in these surveys, and eventually with the Square Kilometre 
Arra>Q 
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